Breakdown of the Mott insulator: 
Exact solution of an asymmetric Hubbard model 



00 

o 



o 



> 

m 
O 

o 

00 
On 
-I— > 



I 

O 

o 



X 
S3 



Takahiro FukuicI 

Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611, Japan 
Norio Kawakami 

Department of Applied Physics, Osaka University, Suita, Osaka 565, Japan 

(December 1, 1997) 



The breakdown of the Mott insulator is studied when the 
dissipative tunneling into the environment is introduced to 
the system. By exactly solving the one-dimensional asym- 
metric Hubbard model, we show how such a breakdown of 
the Mott insulator occurs. As the effect of the tunneling is 
increased, the Hubbard gap is monotonically decreased and 
finally disappears, resulting in the insulator-metal transition. 
We discuss the origin of this quantum phase transition in com- 
parison with other non-Hermitian systems recently studied. 

PACS: 71.27.-ha, 71.30.-)-h, 05.30.-d 



I. INTRODUCTION 

Strongly correlated electron systems in low dimensions 
have attracted renewed interest since the discovery of the 
high-Tc superconductivity 0,^. In particular, systematic 
study of electron systems in the Mott insulating phase is 
believed to provide an important key step to understand 
such unconventional phenomena p\. To clarify correla- 
tion effects on the Mott insulator it is efhcient to 
probe the response of the system to external fields. For 
semiconductors with weak correlations, the charge gap 
due to the one-body band structure is collapsed by strong 
external electric fields, which is well-known as the Zener 
breakdown [Q. For the Mott insulator, however, it is not 
straightforward to study such a problem theoretically, be- 
cause it is essentially a non-equilibrium many-body prob- 
lem. To observe the breakdown of the Mott insulator may 
indeed provide new paradigm for strongly correlated elec- 
tron systems in non-equilibrium conditions. Although 
there have not been many works in this context so far, 
such attempts have been recently started experimentally 

In this paper, we address the problem how such break- 
down phenomena occur for the Mott insulators in non- 
equilibrium conditions. Since it seems still difficult to 
deal with the Zener breakdown directly, we here propose 
a related simple model, and study its properties exactly. 
Namely, as a typical example of the Mott insulator, we in- 
vestigate the one-dimensional (ID) Hubbard model with 
the dissipative tunneling into the environment, which has 
a non-Hermitian character. By using the exact solution 



of the model, we show that there exits a critical strength 
of the tunneling, at which the Hubbard gap is closed and 
the insulator-metal transition occurs. We determine the 
phase diagram by calculating the critical line exactly. 

One of the characteristic aspects in our model is that 
the effect of the asymmetric tunneling gives rise to imag- 
inary gauge potential, which may be regarded as a gener- 
alization of a conventional (real) gauge potential induced 
by twisted boundary conditions |^-|l3|. In this connec- 
tion, we wish to note that such unconventional systems 
with imaginary gauge potential have attracted much cur- 
rent interest for ID random systems [p^| " [2C[| and spin sys- 
tems ||2T|-p9[|. In particular, Hatano and Nelson [Q have 
recently shown that the delocalization transition for the 
pinned flux state of superconductors in cylinder geome- 
try can be described by a ID tight-binding model with 
asymmetric hoppings. They have demonstrated the im- 
portance of the imaginary gauge potential due to the non- 
Hermitian hopping term, which can indeed delocalize the 
"insulating" state caused by the randomness. This type 
of delocalization phenomenon indeed happens even for 
the Mott insulator, as discussed in the next section, and 
both phenomena have a close relationship to each other, 
except that the insulating phase in the present case is 
formed by the electron correlation, not by the random- 
ness. This mechanism is also applied to the collapse of 
the Ising-type spin gap for the asymmetric XXZ model, 
which has been found in ref. 25|. 

This paper is organized as follows. In Sec. II, we intro- 
duce the model Hamiltonian, and mention the relation- 
ship to ID disordered systems with an imaginary vector 
potential and also to the asymmetric XXZ model. In sec. 
Ill, we solve the model by means of the Bethe ansatz 
method, and show how the insulator-metal transition oc- 
curs when the tunneling amplitude is increased. We dis- 
cuss characteristic properties in the insulating phase and 
the metallic phase separately. In Sec. IV, we summarize 
our results. 



II. MODEL HAMILTONIAN 

Let us start with the ordinary ID Hubbard model with 
the on-site repulsive Coulomb interaction U. When an 
external bias voltage is applied between the first site and 
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the last site, it induces non-equilibrium currents, which 
may be essential for the dissipation into the environment. 
In order to effectively formulate the tunneling of electron 
currents into the environment, we introduce asymmetric 
tunneling (or hopping) terms at two end points (1st and 
Lth sites), namely the asymmetry between T]^ic\^^cia 

and TiLc\^CLcr with real amplitudes Til 7^ T^i- For ex- 
ample, when Til < Til, we can effectively describe the 
situation where electrons may be supplied at the first 
site, and may dissipate at another boundary. It may 
be thus expected that this type of asymmetric tunnel- 
ing at boundaries may describe some essential features, 
especially a "direction" for the system, caused by non- 
equilibrium currents. 

In order for the tunneling rate into the environment 
to be a macroscopic quantity, Tli {Til) should be ex- 
tensive. To obtain physically sensible result, it is con- 
venient to parameterize them as boundary conditions, 
CL+ia — e~^*cio- and c|^_|.ig. = e^*c|^, where ^' (order of 
unity) controls the strength of the asymmetric tunneling. 
In this case, we have Tli = — te~^* and Til = — te^*. 
This expression implies that the present boundary con- 
dition can be regarded as a twisted boundary condition 
with an imaginary twist angle. Recall that the ordinary 
twisted boundary condition is equivalent to introducing 
a magnetic flux threaded in the ring system [^,0 . Sim- 
ilarly, the imaginary twisted boundary condition in our 
model effectively induces the imaginary gauge potential 
in a ring system p^ ] . Namely, if we use a gauge transfor- 
mation, the Hamiltonian is now cast into the form with 
periodic boundary conditions. 



j = l cr 
L 



(1) 



where U > 0. It is clearly seen that each hopping term 
acquires imaginary gauge potential The eigenval- 
ues of the Hamiltonian always appear in complex con- 
jugate pairs. Since ^^({-0^}) = H{—{-ipj}), the right- 
eigenvector of H{{ipj}) is equivalent to left-eigenvector 
of H{—{ipj}) with the same energy. 

It is to be noticed that the Hamiltonian is non- 
Hermitian, because we have eliminated the environment 
by introducing the asymmetric tunneling instead, which 
as a result causes dissipative effects. The non-Hermitian 
property may be natural since we are now concerned with 
non-equilibrium phenomena. In the gauge-transformed 
representation (0), the hopping term feels an imaginary 
gauge field ^ at every site. Now the relationship to 
the ID disordered tight-binding models with an imagi- 
nary vector potential |p^-p^ and the asymmetric XXZ 
model [pT[-p9t is clear: The hopping term in Eq. (|l|) 
takes essentially the same form as found in such models. 



The difference among them is the mechanism to realize 
the insulating phases when = 0: The Mott insulator 
is stabilized by the Hubbard interaction for the present 
model, whereas the Anderson insulator is realized by ran- 
dom potentials for disordered tight-binding models, and 
the spin-gapped state appears by the Ising-anisotropy 
for the asymmetric XXZ model. Therefore, the present 
insulator-metal transition may have essentially the same 
origin as for the Hatano-Nelson model |Q and also the 
asymmetric XXZ model p^ , p8| |. 

In the next section, we show that the charge gap of the 
present system closes due to the asymmetric tunneling 
effects, and consequently the Mott insulator breaks down. 



III. BETHE ANSATZ EQUATIONS 

Let us now turn to the solution of the above Hamilto- 
nian. The Bcthc ansatz method still works even for such 
unconventional Hamiltonian, because asymmetric hop- 
pings in the model are incorporated as twisted bound- 
ary conditions with imaginary twist angle ^1 32 1, as 
mentioned above. In fact, the asymmetric XXZ model 



have been extensively investigated by the Bethe ansatz 
method |2l|-[29|]. It has been shown that the gap due 
to Ising-anisotropy closes when the imaginary twist is 
increased Now, applying the Bethe-ansatz tech- 

nique to our model, we diagonalizc the Hamiltonian (|l|) 
to obtain the Bethe-ansatz equations |10 3^-|35[|, 



kjL = 2TrIj + iL^ + ^ 6* (sin - A^), 

f}=i 

^ 0(sin k, -K)^2nJ^-J29( ) , (2) 



0=1 



where the two-body phase shift is 0{x) = — 2 arctan(x/u) 
with u = U/{At). Here the quantum numbers are 
I J = Ni/2 (mod 1) and J„ = {N - + l)/2 (mod 
1). The branch cut of 9 is set along (— ioo, — ic) and 
(+ic, +ioo). The energy of the whole system is given by 

E — -2t Y^^=i kj . 

In what follows, we consider the half-filling case, where 
the Mott insulating phase realizes for any finite U when 
= 0. Hence, we set the number of electrons N — L 
and the number of down-spins — L/2. In Fig. |^, we 
show the distribution of the rapidities on the complex k 
and A plain, which are calculated directly by Eq. (||) for a 
small system. At ^ = 0, all rapidities are on the real axis, 
whereas they form nontrivial curves on the complex plane 
for a finite 4". Their distributions are always symmetric 
with respect to the imaginary axis, which ensures that the 
ground state energy is real. Though imaginary twisted 
boundary conditions are imposed on the charge sector, 
the spin distribution is also affected slightly. 
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FIG. 1. Ground-state distribution of rapidities for the 
1/ = 38 system (half filling) with u = 1 and 
h = 0.8 (* = 0.242797). The system with u = 1 has 
b^^ = 0.881374 (*cr = 0.246349). The solid curve is the 
distribution of charge rapidities in the bulk limit calculated 
numerically by Eq. ^ . The horizontal and vertical axes are, 
respectively, real and imaginary part of k and A. 



A. Insulating phase 

In order to study critical properties of the model, we 
should carefully observe the behavior of rapidities. To 
this end, let us define the counting functions, as usual, 
Zc{k) and Zs{X) satisfying Zc{kj) = Ij/L and Zs{Xa) — 
Ja/L, respectively. In the bulk limit, they are given by 



Zc(fc) = 



k 
1 

1 



1 



dke{smk~ X)p{k) 
A - A' 



dAe'(sinfc- A)cr(A), (3) 



a(A'), 



(4) 



where the distribution functions arc p{k) = z'^{k) and 
(7(A) = z'^{X), satisfying 



p{k) = 



1 

2^ 



cos k 

271 



dXe'isink- X)a{X), 



ct(A) = -^ / dke' {sin k - X)p{k) 
27r Jc 



A/^.A'.'(A_^).(A'), 



(5) 



(6) 



and where C and S denote the curves on which the solu- 
tions of the rapidities lie. 

Now, starting from the Mott insulating phase realized 
for small we wish to deduce the critical point at which 
the Mott insulator breaks down. For this purpose we 
need to study the analytic properties of the ground state. 
First, we denote the end points of the curve C as a + ib 
with a = ±7r in the bulk limit, whereas those of S are 
±cxD. In the complex k plane, the poles of 6''(sinfc — A) 
in Eq. (H) nearest to the real axis are ztiarcsinh u (mod 
tt). Therefore, if & < 6cr, where 



C and S can be deformed, respectively, as — tt + ib ^ 
— TT —^^T—^^^ + ib and ~oo — > oo with straight segments. 
Then the distribution functions are given by the analytic 
continuation of those with 5* = 0, which means that the 
system is still in the insulating phase. For reference, we 
plot the curve C calculated by Eq. (^) in the bulk limit 
as a solid line in Fig. 0, from which one can observe that 
rapidities for a finite-size system sit exactly on C of the 
bulk limit except for the end ones which deviate slightly. 

Since b characterizes the curve C, there is one-to-one 
correspondence between ^E* and b. By using Zc{k = ±7r-|- 
ib) = ±1/2 in Eq. (^), which is valid by definition, we 
obtain the relation between ^ and b, 



*(&) 



dA6l(A-|-isinh6)CT(A). 



(8) 



Eventually, we can determine the critical value of ^! by 
the formula 'i'ci ~ hmb_+f,cr-o ^(&) with Eq. (Q), at 
which the analytic property of the ground state changes, 
namely, the breakdown of the Mott insulator occurs. In 
Fig. H, we plot the critical values of ^I^cr as a function of 
u. For small u, the Hubbard gap is exponentially small, 
and hence it is easily collapsed by small 5*. For large u, 
we have the dependence of ^'c arcsinh m, so that we 
observe the crossover behavior around u = I. 
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FIG. 2. Phase diagram for the insulator-metal transition. 
The solid line is the critical line ^ = ^'cr(ii). 



1. Hubbard Gap 

Now that we have determined the critical ^'c in terms 
of singular properties in the ground state, let us di- 
rectly evaluate the Hubbard gap to confirm it. At 
^ = 0, the charge excitation has the Hubbard gap 
for any finite U . We will show that the gap becomes 
smaller as is increased, and is finally collapsed at 
^! = The Hubbard gap is here defined Q by 

A = [/ — 2/i_, where the chemical potential p,^ is given 
by p- = E{N'f,Ni) - EiNf - l,iV|). To obtain /i_ in 
terms of a one-hole excitation energy, let us make a hole 
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state at Ih in consecutive quantum numbers Ij . Then we 
readily find that the energy change e(fc/i) = E — Eg,s. is 
given by the analytic continuation of that with = 0, 



e(fc,0 = 2t 



COSfc/i+ / dcj -; „ I I, 



(9) 



where corresponds to the rapidity of the hole Ih ■ The 
lowest charge excitation at kh — ±7r + ib determines the 
Hubbard gap. Therefore, setting /i_ = — e(±7r + ib), we 
reach the expression for the Hubbard gap for finite 



A(6) = 4t 



cosh b - 



; sinh b 



duj- 



I) 



(10) 



where Jn{(^) is the nth Bessel function. It is now checked 
that the Hubbard gap is a monotonically decreasing func- 
tion of which confirms that the dissipative tunnelings 
into the environment indeed deduce the Hubbard gap. To 



evaluate the gap at 5* 
the expression, 



it is convenient to rewrite 



n=l 



A(6)/4t = u-cosh6-^(-)" Vl + (2ufi- sinh 6)2 

(11) 



+ + {2un + sinh 6)2 — Aun 



It is readily seen that A(6cr) = 0, by setting sinh6cr = u. 
Namely, at ^E" = ^'cr which has been determined by the 
analytic property of the ground state, the Hubbard gap 
indeed closes. Note that the expressions for the gap in 
this subsection is vahd only for < ^fcr, i.e., 6 < 6cr. 
Therefore, we can confirm that on the critical line in 
the phase diagram Fig. ^ the breakdown of the Matt 
insulator occurs due to the effect of dissipative tunnelings. 



2. Ctiarge excitations 

Next let us discuss how densely the charge excita- 
tion is distributed above the Hubbard gap. It is usu- 
ally discussed by introducing the distribution function 
for charge excitations, D{e) — dzc/de, where e{kii) = 
4t[c + e{kh)/2t] with Eq. (|), taking the gap into ac- 
count. For the present model, however, it may be more 
suitable to define it as 



Die) 



Rep 
iRee'l 



(12) 



since the energy of the elementary hole excitation is com- 
plex in general. Here, we evaluate p and s along the curve 
C in the bulk fimit calculated by using Eq. (^), e.g., the 
solid-curve in Fig. 0. We shall refer to -D(e) as the density 
of states (DOS) for simplicity. 

In Fig. ||, we show DOS ( p^ ) for u = 1 case. One can 
see that at the band edge, DOS diverges as is the case 
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FIG. 3. Density of states (solid-line) for charge excitations 
and imaginary part of the excitation energy (dashed-line) for 
u — 1 system with b=0.1 (a) and b=0.8 (b). The correspond- 
ing * is 0.0399396 for (a) and 0.242797 for (b), respectively. 
We here set 4t = 1. The straight line is a guide for gaps. 



for the conventional ID Hubbard model. For reference, 
we also plot the imaginary part of the excitation energy 
e. Here, since Im£(— fc/i) = —Ime{kh), we plot —lnie{kh) 
for kh > 0. For a small [Fig. ^ (a)], eigenvalues are 
almost real, and behavior of the DOS is quite similar to 
^I' = system, as should be expected. As ^I' increases the 
charge gap certainly decreases. For a large 4" near "ifcr 
[Fig. I (b)], DOS still diverges. 

Now we observe the properties around the band edge 
in a bit detail. To this end, we need to calculate p and 
e' at fc = TT + z6cr- In a similar way to Eq. (pTf), we have 



27rp(7r + ib)^l~ cosh6^(-)"+i 

n=l 

1 1 



■yj{2un - sinh 6)2 + 1 ^(2un + sinh 6)2 + 1 
Therefore, 

p{ti + i6cr) = 0, 

p'(7r + i6c,) -^(li2 + l)F(w)/^, 



(13) 



p"(7r + ib„) ^ Zc^Ju^ + 1F{u)/tt, 



where 



F{u)^j2(~y' 



{2n + l)u 



^(2n+ 1)2^2 + 1 



(14) 



(15) 



Here we calculate p' and p" for later convenience. We 
now find that near the critical point p is real and ap- 
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proaches to as p{n + ib) —ip'{-K + i6cr)(^cr — 
Therefore, at first sight, the DOS defined by ( |l^ ) does 
not seem to diverge at the band edge. However, noting 
e{-K + ib) = A (6), we have 



e'{TT + i6cr)/4< = iu~ iy/u'^ + 1 



div- 



1 



r)2u\uj\ 



(16) 



Namely, e'(7r + ibcr) is pure imaginary. Therefore, the 
DOS of the system near the critical point diverges at the 
band edge, as is seen in Fig. ^ (b) [so far as we use the 
definition (|l2|)], although p approaches 0. Anyway, the 
excitation above the Hubbard gap is unconventional with 
complex energy, which reflects the fact that the present 
system possesses non-Hermitian tunneling to the envi- 
ronment. Thus, when the Mott insulator breaks down, 
the resulting metallic state would have a dissipative char- 
acter with complex eigen energy. 



B. Metallic phase 

So far we have shown that the Mott insulator breaks 
down at ^' = ^cr- Now the question is what kind of state 
is realized beyond the critical "ifc- We shall show that it 
is a dissipative metallic phase without charge excitation 
gap. Let us recall first that in the Mott insulating phase, 
the real part of rapidities completely fills up the lower 
Hubbard band in the region [— tt, tt]. In order for the sys- 
tem to enter the metallic phase, the "Fermi points" ±a 
for the real part of the rapidity should leave away from 
±7r, when we increase 4* > 'Sc. Indeed, from the expres- 
sion (p^ ) we see that p{tt + ib) decreases if we increase 
b and eventually at & = 6cr, p{t^ + *^cr) — 0. Hence we 
have a < TT when b > 6cr- To explicitly determine the 
Fermi points a (and b as well), we deform the curves C 
and S as the straight segments —a + ib a + ib and 
— cx) -|- iXi — > oo + iXi , respectively. Consistency between 
Eq. (^ and JgdXa{X) = 1/2 requires 



sinh b — u < Xi < cos a sinh b + u. 



(17) 



Similarly, from Eq. (|^) and dkp{k) = 1 [ or from Eq. 
(0) ], we have. 



a 1 
TT 27r 



dXa{X)\e{iiin{a + ib) - X) 
(sin(-a-f i6) - A)|. (18) 



By using Eq. corresponding ^ is given by 

/OQ-\-iXi 
dXe (sin(a + ib) - A) a{X). 
-oo+iAi 

(19) 

These are the basic equations to determine the effective 
"Fermi points" ±a. 



Let us solve these equations explicitly for two limiting 
cases. First, consider the large 4* region, namely the 
limiting case of strong dissipation. When b — )■ oo, we 
have a — > with < a^e''/2 < 2u from the relation (|lT 
and hence sin(±a + ib) ^ X — iXi ^ zLae'' — A -I- iu' 



±oo 



where —u < u' < u due to Eq. (fl 



Then 



e (sin(±a + ib)-X- iXi) T[7r + 2u/ {a^e'')]. Finally, 
from Eqs. (^8|) and (p^), we end up with 



u 

TT 



^a^e"^-, VI/ = 6. 



(20) 



This equation states that for a large vl>, a indeed ap- 
proaches in the form a cx e^*/^. Therefore, we can say 
that a large enough vp effectively drives the system to the 
low-density limit of the ordinary Hubbard model except 
that the energy eigenvalues are now complex reflecting a 
dissipative metallic state. 

Next consider the metallic phase just beyond the crit- 
ical point. Setting a — tt ^ 5a and 6 = 6cr + Sb, we 
have 2VvFTTSb < u{da)'^/2 from Eq. (|l^. Hence 
we see 6b = 0{6a^). Similarly, setting p{k + ib) — 
Pcr{k + ib) + Sp{k) and cr(A + iXi) = crcr(A + iAi) + Sa{X), 
where pa- and dcr are distribution functions on the critical 
point, we immediately find 6a{X) — 0{6a^) by expanding 
Eqs. (^) and (^) in power series of Sa and Sb. Then Eqs. 
(|l|) and dH) reduce to 



da — dXacr{X) 
(5* = i5a + 5b 

+i / (iA(Tcr(A) 



'6a + ie"6a6b + -6"' 6a'^ 
6 



0(^4), 



-e'6a + -jQ"6a^ 



+ 0((Sa3), (21) 



where 6^ = ^ - and d' = d6{sink - 

Tr+i{b„-o)- Contributions from 6a are included 
in higher order terms. A little algebra leads us to the 
final relations 



6b ^ 



6hiPcr(7I" + «fecr)] 



6a^ + 0{6a^), 



6^ = <l+'-^^^^^^±^]6b + 0{6a% 



(22) 



where p' and p" are given by Eq. (pT[). 

The flow of the end points is schematically illustrated 
in Fig. ^ as a function of vp. We can now see how the 
Hubbard gap disappears at half-filling: The larger vp is, 
the longer the curve C becomes (see Fig. |^). As a result, 
the density of charge rapidities at the end points (±7r) 
decreases and vanishes at = vj/^^. [gee Eq. (|l^]. If 
we further increase ^, the end points (Fermi points) are 
forced to be away from tt, realizing the metallic state. It 
is to be noted again that this metallic phase has a dissi- 
pative character with complex eigenvalues, which shows 
sharp contrast to the ordinary metallic phase. 
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FIG. 4. Schematic illustration of the behavior of the end 
points ±a + ib. They move toward the direction indicated by 
arrows if we increase ^. 



IV. SUMMARY 

In summary, we have proposed a non-Hermitian Hub- 
bard model to study the insulator-metal transition, which 
may describe a many-body analogue of the Zener break- 
down, driven by asymmetric tunneling into the environ- 
ment. By solving the model exactly, we have shown 
that such a tunneling indeed causes the breakdown of 
the Mott insulator, and drives the system to a metallic 
state. We have studied the basic properties of the in- 
sulating phase and the metallic phase, and clarified the 
mechanism of the transition. In the sense that this tran- 
sition is caused by the imaginary gauge potential induced 
by the asymmetric tunneling, it shares essential proper- 
ties with other interesting non-Hermitian systems such 
as the model of Hatano and Nelson for the delocaliza- 
tion of pinned flux state in superconductors and also the 
collapse of the Ising-type spin gap for the asymmetric 
XXZ spin chain. Although the experimental study on 
the breakdown of the Mott insulator has not been done 
systematically so far, we think that this may provide an 
interesting subject in the correlated electron systems in 
non-equilibrium conditions. 
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